Persistent currents in a Bose-Einstein condensate in the presence of disorder 
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We examine bosonic atoms that are confined in a toroidal, quasi-one-dimensional trap, subjected 
to a random potential. The resulting inhomogeneous atomic density is smoothened for sufficiently 
strong, repulsive interatomic interactions. Statistical analysis of our simulations show that the gas 
supports persistent currents, which become more fragile due to the disorder. 
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I. INTRODUCTION 

Superfluidity is one of the most fascinating problems in 
physics. Superfluidity includes a whole class of phenom- 
ena that appear in various quantum systems such as 
liquid He, superconductors, and nuclei. More recently, a 
lot of effort has been put on the study of superfluid prop- 
erties of dilute vapors of trapped atoms, including, for ex- 
ample, vortex states @, H, 0,111], transport properties @, 
etc. Although these gases are very dilute, they are super- 
fluid because of their extremely low temperature. Their 
diluteness makes cold gases of atoms an ideal system for 
testing superfluid properties. For example, their theoret- 
ical description is much easier as compared to other su- 
perfluids, like e.g., liquid Helium. Furthermore, they can 
be manipulated in numerous ways, including the form of 
the trapping potential, and their coupling constant. 

In the present study we examine the stability of persis- 
tent currents in trapped superfluid atoms that are con- 
fined in toroidal, quasi-one-dimensional traps. While 
conservation laws imply trivially the conservation of mo- 
mentum/angular momentum even in a moving/rotating 
classical system, the crucial question is the robustness 
of the current-carrying state(s) against impurities or 
anisotropies; it is exactly this feature that distinguishes 
a quantum from a classical system. 

Many theoretical studies have examined this problem, 
see, e.g., Refs. 0, 1, 1 M, EH, 03 03 • Experimentally, 
numerous recent studies have considered this question in 
connection with trapped atomic gases [H, 03 03 Ell- 
in another class of experiments - which are closely con- 
nected with the present problem - the propagation of 
Bose-Einstein condensed atoms in mag netic waveguides 
has also been examined thoroughly [19(. 

Interestingly, the problem that we study here is not 
only a toy model that provides answers on some non 
trivial questions, but also it is directly applicable to 
experiments that have created Bose-Einstein conden- 
sates of atoms in toroidal traps [2(j. Furthermore, 
the stability of supercurrents in these mesoscopic sys- 
tems may have revolutionary technological applications 
[2l| : similar arguments may also be applicable to other 
(nanoscale/ mesoscopic) systems, such as quantum rings 
and quantum dots [22j. 

In what follows we first describe our model in Sec. II. 



In Sec. Ill we consider the effect of a random potential 
on the atomic density of a static cloud. In Sec. IV we 
discuss how one may understand the existence of per- 
sistent currents as a result of the metastability of the 
current-carrying state. Section V presents the results for 
the effect of disorder on the stability of superflow within 
the mean-field approximation. In this section, we con- 
sider many random potentials and analyze statistically 
our results. In Sec. VI we solve the same problem within 
the Bogoliubov approximation. Finally, in Sec. VII we 
discuss our results and give some general conclusions. 



II. MODEL 

In our model we consider a tight toroidal trap, which 
freezes the transverse degrees of freedom of the motion of 
the atoms. Further, we assume for the atom-atom colli- 
sions the usual contact potential, Vint(r— r') = Uq5(t~ r') 
with Uq = Air?i 2 a sc /M . Here, a sc is the scattering length 
for elastic atom-atom collisions and M is the atomic 
mass. Since the transverse degrees of freedom are frozen, 
the three-dimensional field operator may be written as 
t/N/(RS) 9(6), where 9(6) is the field operator associ- 
ated with the motion of the atoms along the torus. Here 
9 is the azimuthal angle, N 3> 1 is the atom number, R 
is the radius of the torus, and S is the cross section of 
the torus (with R \/~S). Therefore, the Hamiltonian 
of our quasi-one-dimensional system is (23l . l2~ij 



H/N = / ¥(9) 



d 2 



2MR 2 d9 2 



V(9) 



9(6) d9 



^nnoUo { ¥(6)¥(6)9(6)9(6)d6, (1) 



where V(9) is the external potential, which is taken to 
be piecewise constant and random (see bottom graph in 
Fig. 1). Finally, no is the average atom density, uq — 
N/(2nRS). 

Within the mean-field approximation the condensate 
order parameter 9 satisfies the nonlinear equation 



d 2 9 



V(6)9(6) + 2TT-f\9(9)\ 2 9 = ii9, 



(2) 



d9 2 

where we have set Ti = 2M = R = 1; V(ff) and the chemi- 
cal potential [i are measured in units of the kinetic energy 
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T = h 2 /(2MR 2 ). The ratio between the interaction en- 
ergy and the kinetic energy is equal to 7 = n Uo/T. This 
parameter is also equal to 7 = ANa sc R/ S . 

The mean-field approximation is valid when the di- 
mensionless quantity 7 is -C N 2 [25]. Further, the limit 
of weak interactions is achieved when 7 < 1, while the 
Thomas- Fermi limit is achieved when 1«7< iV 2 . For 
the typical values of 7 that we consider in the present 
study, the mean-field approximation is applicable, and 
the healing length £ is on the order of the radius of the 
torus, since £/R = 7~ 1//2 . The opposite Tonks-Girardeau 
limit of impenetrable bosons is achieved when 7 becomes 
of order N 2 [2j|, where £/R ~ 1/N. 



III. DENSITY PROFILE IN THE PRESENCE 
OF A RANDOM POTENTIAL 

As a first step, we consider the static profile of the gas 
in the presence of the random potential V(0). We choose 
the length scale d of variation of V(8) to be R/10, as in 
the experiment of Ref. [15j. In this experiment, the axial 
size of the condensate was « 110 //m, the radial size was 
w 11 /im, the smallest length scale of each "speckle" was 
« 10 //m, and the average distance between speckles was 
« 20 jum. 

There are four energy scales in our problem, namely 
the kinetic energy associated with the torus % /(2MR 2 ), 
the characteristic depth of each separate potential Vb, 
the zero-point kinetic energy Ti 2 /(2Md 2 ), and the typi- 
cal interaction energy UqUo- Since we choose d = R/10, 
therefore H 2 /(2Md 2 ) is two orders of magnitude larger 
than h 2 /{2MR 2 ). We also choose V to be of the same 
order as K 2 / (2M R 2 ) , which implies that the lowest state, 
as well as the low-lying excited (eigen)states of the non- 
interacting problem are localized around the minima of 
the random potential, but there is also a significant over- 
lap between the maxima in the density, as shown in Fig. 1 
for 7 = 0. If one chooses the zero of the energy to be at 
the maximum of the random potential, quantum mechan- 
ics implies that in the effectively one-dimensional prob- 
lem that we consider here there is always at least one 
bound state, and thus at least one exponentially local- 
ized/delocalized state [26l[27j. 

Under these conditions, depending on 7, there are 
three different regimes. For zero/ weak interactions, 7 <C 
1, the order parameter is equal/close to the lowest eigen- 
state of the external potential. As 7 increases, the den- 
sity variations get suppressed. When 7 > 1, the density 
becomes homogeneous, as shown in Fig. 1. 



IV. STABILITY OF PERSISTENT CURRENTS 

We turn now to the stability of persistent currents. 
Before we consider any disorder, it is instructive to see 
how one understands the existence of persistence currents 
for a constant potential V(6). 
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FIG. 1: The density j*(6»)| 2 of Eq. © for three values of 
the coupling constant 7 = 0.0, 1.0, and 10.0, for the random 
potential V{&) shown in the bottom graph. 



The basic idea is that for sufficiently strong interac- 
tions, the dispersion relation £{l), i.e., the energy per 
atom as function of the angular momentum per atom, 
develo ps a barrier between the states with I = 1 and 
I = |28l . I29I l30l ]. Recalling that when 1 = 1, there is 
a vortex state that is located at the center of the torus, 
while for I = the vortex is at an infinite distance away, 
in order for the vortex to exit the torus, there has to be 
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FIG. 2: The dots show the average value (7c) and the bars 
show the standard deviation <x(7 c ), of the critical coupling 
7c versus the "strength" of the random potential k. These 
values are calculated from the examination of 1000 different 
random potentials. For each specific potential, the values 
of V(8) at each subinterval are chosen independently from a 
uniform random distribution. The solid curve shows 7 C for the 
specific potential shown in Fig. 1. The dashed curve shows 7 C 
that results from Eq. (j4|, i.e., from the toy model described 
in the text, for the same potential of Fig. 1. 



a node in the density of the atoms. However, this node 
costs interaction energy, if the coupling is strong enough 
(since for repulsive interactions the interaction energy is 
minimized for a homogeneous density). The energy bar- 
rier that separates the rotating state with I = 1 from the 
lowest-energy state with I = thus allows for the exis- 
tence of persistent currents, making the timescale for the 
decay of the current exponentially long. 



V. PERSISTENT CURRENTS WITHIN THE 
MEAN-FIELD APPROXIMATION 

We start our analysis of the stability of persistent cur- 
rents within the mean-field approximation, for relatively 
strong disorder, where Vq ~ h 2 /(2MR 2 ) ~ tiqUq. In 
the presence of any spatially-dependent potential V(9), 
the angular momentum is not a good quantum number, 
and even more so in this case, where the disorder is not 
treated perturbatively. Still, one may examine the ener- 



getic stability of a state with one unit of circulation 31 1 ; 



energetic stability also guarantees dynamic stability 
Clearly such a state can only be metastable and not the 
lowest-energy state of the system, since the non-rotating 
state has a lower energy. 

Before we describe our numerical results, it is instruc- 
tive to present a toy model that describes the problem 
qualitatively. We thus consider the random potential 
shown in the lowest graph of Fig. 1, multiplied by some 
dimensionless constant K, i.e., nV(6); essentially k is the 
"strength" of the disorder. Although this argument can 
be generalized, for the sake of simplicity we consider a 
(very limited) truncated order parameter 



that has an expectation value of the angular momentum 
per atom equal to I. Here $ m = e" n6 /V2n and A is a 
phase that is determined below. This order parameter 
gives an energy per atom that satisfies the equation 



£_ 

N 



7 



kV = (1 + 7)/ - ~fl 2 - 2K\V e \y/l(l - I), (4) 



(3) 



where V c = f ^ V {9) cos 9 d9 / (2tt) , and V = 
J* V(9)d9/(2n) is the zero component of the Fourier 
transform of V(9). If V c is negative/positive, the phase 
A is chosen to be A = 0/V, in order for the last term in 
Eq. (0| to be always negative. Since, according to Eq. Q, 
when k = 0, £ (I) develops a local minimum at I = 1 for a 
value of 7 = 1, the critical value of the coupling 7 C that 
gives metastability of the current is 7 C = 1 for the spe- 
cific ^»tr- According to Eq. for K ^ 0, the last term 
may destabilize the current, forcing 7 C to increase with 

as shown in the dashed curve of Fig. 2. In other words, 
the coupling has to be sufficiently strong in order for the 
state with nonzero circulation to be stable. Clearly this 
curve is qualitatively, but not quantitatively correct. 

Turning to the numerical results, we attack this prob- 
lem via the method of imaginary time propagation [33| . 
We choose different realizations of random potentials of 
the form shown in Fig. 1, multiplied by n. For each spe- 
cific potential, the values of V(9) at each subinterval are 
chosen independently from a uniform random distribu- 
tion. For each specific V(9), we identify the critical value 
of the coupling 7 C that is necessary to obtain stability for 
a state with unit circulation for a specific k, starting with 
$1 as an initial condition. The solid curve in Fig. 2 shows 
7c(k) for the specific random potential shown in Fig. 1. 
For k — ► 0^ , our results imply that j c — 3/2; this result 
is analyzed further below. We also perform a statistical 
analysis of the values of 7 C that we get using this method, 
i.e., we compute the average value ("y c ) and the standard 
deviation <r(7 c ), for 1000 different random potentials for 
each value of k. Figure 2 shows the result of these cal- 
culations. The values of 7 C that we get from the random 
potentials that we consider are approximately Gaussian 
distributed. 

The states that come out of our calculation are differ- 
ent than $1 or $0, as they are not homogeneous, because 
of the (random) potential. Figure 3 shows the calculated 
density of the gas, for a circulation equal to zero (solid 
curve) and unity (dashed curve), for the specific random 
potential shown in Fig. 1 (with k = 1). To get these 
densities, we use the initial condition $0 = 1/ \[2/k and 
$1 = e i0 /\/27r, respectively. 



VI. BOGOLIUBOV APPROACH FOR WEAK 
DISORDER 

Another way to attack this problem, incorporating the 
interactions, as well as the disorder, is via the Bogoli- 
ubov transformation (24|. This method gives the whole 
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FIG. 3: The density |*(0)| 2 of the cloud that results from the 
imaginary-time propagation, for 7 = 5, for the specific poten- 
tial of Fig. 1, with k = 1. The solid/dashed curve corresponds 
to the state with zero units/one unit of circulation. 



excitation spectrum and also the depletion of the con- 
densate, which is assumed to reside at the single-particle 
state $1 = e i6 /y2w. This assumption, as well as the as- 
sumption for the existence of a Bose- Einstein condensate 
requires that the disorder is not too strong. 

Following the usual tricks, we replace the creation and 
annihilation operators of atoms with angular momentum 
m = 1 , c\ and c\ , with \f~N \ , where N\ is the number 
of atoms occupying the m — 1 state. Then, since TV = 



j, one can write for the Hamiltonian 



H — N — -/N/2 =^(m 2 -2m + l)c\_ m ci- 



m>0 



+ (m 2 + 2m + j)c\ +m c 1+m + 7(ci_ m ci+ 



"'' "T" c I-m C I+m) 



+e* m (cx- m + 4+m) + e m(c\- m + Cl+ m ),(5) 



In order for the excitation spectrum to be positive, 
the coefficient of the operator a^am implies that 7 > 
2 — m 2 /2. The most unstable mode is the one with m = 
1, and therefore the critical coupling constant for the 
existence of persistent currents is 7 C = 3/2. This result 
was first derived in Refs. [23|, (24|, within the truncated 
basis set of the states with m = 0, 1, and 2. According 
to the argument given above, this is the exact value of 
7 C for the stability of persistent currents, as we have also 
found numerically within the mean-field approximation. 

The eigenstates \n a , m , np tT n) of the diagonalized 
Hamiltonian of Eq. ([7]) (i.e., of the number operators 
a^am and fi^fim) carry an angular momentum L = 

n-y: 



m>0 



m(n a ^ m - n/3, m ), i.e., L — N ± 0(1). From 
Eq. ([7]) we also see that the disorder potential enters the 
Hamiltonian only via the last term, which can also be 
derived perturbatively. As a result, one cannot see the 
dependence of 7 C on k within this method. In order to see 
this dependence, one would have to consider a conden- 
sate of non-uniform density, or equivalently a condensate 
with at least two single-particle states macroscopically- 
occupied, e.g., of the form 10^°, l Nl ), as in Eq. ©■ 

Equation ([7J can also give the depletion of the 
condensate J2 m >o( n a,m = Q,n^ m = 0|c}_ m ci_ m + 
c[+m c ^+m\ n a,m = 0, n^ jm — 0), which consists of two 
parts, the one resulting from the interactions and the 
other resulting from the disorder, 



A/V ir 



m>0 



2vi 



7 + m 



disorder 



= 27V 



= V 

m >o m^/m 2 + 2j 

\V m \ 2 (m + 2) 2 



1. 



' m 2 (m 2 + 27 — 4) 2 

m>0 \ ' I J 



(8) 



where e m = y/NV m . The last four terms describe the 
processes where atoms scatter between the condensate 
and the states $i± m because of the (random) potential, 
whose Fourier transform is denoted as V m . We define the 
new operators 

_ | (m + 2)u m + (m — 2)v m 

+ m(m 2 + 27-4) 

(m — 2)u m + (m + 2)v m 



m(m 2 + 27 — 4) 



(G) 



and set u rn = cosh# m and v m = sinh# m , with 
tanh(26> m ) = j/(m 2 + 7), in order to eliminate the off- 
diagonal terms. Then, the Hamiltonian is written in the 
diagonal form 



H — N — 7JV/2 = (~ 2m + m V m2 + 2 7)a I t „a r , 



m>0 



+ (2m + my/m 2 + 2-y) ^ m j3 m 
+m^ 2 T^-(m 2 +7)- A - (7) 

771 +27 — 4 



VII. CONCLUDING REMARKS 

According to the results of the present study, one may 
conclude very generally that disorder destabilizes states 
with nonzero circulation. Physically, any kind of disor- 
der results in an inhomogeneous density distribution of 
the atoms. As compared to the density no of the atoms 
in the absence of any disorder (which is homogeneous for 
repulsive interactions), the minimum and maximum den- 
sity, 77 m j n and 77 max , in the presence of disorder is thus 
both higher, as well as lower than no, n m j n < no < n max , 
as the conservation of the number of atoms implies. 

In other words, in the presence of disorder, there ex- 
ist points along the torus where the local density n(9) is 
lower than 770. These points (where the density is lower 
than no) give the cloud the chance to get rid of its cir- 
culation at a lower energy expense as compared to the 
homogeneous case. Higher values of n enhance the in- 
homogeneities in the density, and thus make it easier for 
the vortex to slip out of the torus. Therefore, the higher 
the strength of the disorder, the higher the interaction 
strength that is necessary to make the compressibility 
high enough, in order for the current to become stable. 
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The above arguments are rather general and apply to any 
external potential. 

The situation we consider here may be realized experi- 
mentally in two different ways: (i) the gas is prepared in 
a rotating state at a temperature above the condensation 
temperature, and then is cooled down to zero tempera- 
ture or (ii) starting from a gas at zero temperature, 
a phase is imprinted 34] . Finally, one could get evidence 
for the presence of current/circulation in such systems 
by, for example, interference experiments (35j . 
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